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Abstract 

We investigate whether the total number of electrons per colli- 
sion, as a function of beam energy, is a potential quark-gluon-plasma 
(QGP) signature. At high beam energies, we find that, after experi- 
mental removal of the noise induced by Dalitz decays, the number of 
produced electrons is increased many times fold from the existence of 
the QGP. The robustness of the potential signal to differing theoretical 
assumptions makes it an attractive experimental parameter. 
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Work is in progress to construct a relativistic heavy-ion collider to test the 
theoretical prediction [1|] that QCD (quantum chromodynamics) undergoes a 
phase transition to unconfined quarks and gluons at high temperatures. Suc- 
cess or failure rides on the identification of an experimental signature which 
indicates when a QGP has been achieved. The complexities of rhic, however, 
are such that a QCD description from first principles will be a long range 
task. Questions of chemical equilibrium, quark effective masses, different 
flavor and gluon temperatures are just some of the many particulars which 
have no universal theoretical consensus. Even a parton description is not 
as adequate for rhic as it is for, say, deeply inelastic scattering, since the 
nuclear medium modifies the parton distribution functions and interaction 
amplitudes. Thus from the experimental point of view, a QGP characteristic 
is desired which is independent of the inclusion of the multitude of indeter- 
minate theoretical fine minutiae. Original ideas to detect the QGP phase 
are based on distributions of produced leptons and photons 0, and to date 
much work is being done to understand the issues involved^]. However, 
what we seek here is a robust QGP signal, that will survive the fine details. 
We demonstrate that the total number of direct electrons per collision as a 
function of beam energy is just such a signature. 

The Result 

In Fig. 1 we represent the simulation of the production of electrons in rhic, 
after the subtraction of Dalitz decays of the low mass hadrons 7r°, 77, p, u, rf, 
(ft, i. e. the directly produced electrons. The Dalitz decays constitute the 
background (noise) from which the signal is to be extracted. While the 
number of mesons expected per collision numbers in the thousands (see Fig. 
2), the accompanying number of directly produced electrons, Fig. 1, can be 
counted by the fingers of one's hands. Identifying and subtracting the Dalitz 
electron pair background does not represent severe difficulties. 

There are two important physical effects associated with the presence of 
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a QGP in conjunction with the measurement of the total number of directly 
produced electrons per collision, as the beam energy is stepped up. When 
the energy density of the hadron phase first reaches the QGP transition 
value, allowing a mixture of the hadronic matter to coexist with the QGP, 
the number of electrons increases more rapidly, with a discernable gradual 
change in slope around the critical beam energy (10-15 GeV/u; see Fig. 1): 
this constitutes the first QGP electron signature. Indeed, if one assumes no 
QGP, from a low beam energy region (say, 5 GeV/u) to a high one (say, 30 
GeV/u), the slope of the lepton yield increases by a factor of 30, while if 
there were a QGP formation at a critical value in between, the slope would 
increase by a factor of 200. This is 6 times larger than in the previous case. 
The fact that there is this change in slope is independent of the theoretical 
uncertainties and is due to the increased degrees of freedom associated with 
the QGP. Past the onset of the appearance of the QGP, as shown in Fig. 1, 
the now more rapid electron production rate results in a value for the number 
of directly produced electrons per collision many times larger than that in a 
pure hadron gas: this is the second QGP signature. Again, physically, the 
QGP produces more direct electrons due to the increased number of degrees 
of freedom. The simplicity in the origin of this signal is at the same time 
its robustness: fine details in the model or uncertainties in the parameters 
used cannot cause a misleading signal. Indeed, we have varied the values 
of different parameters of the model over reasonable ranges and found no 
significant change in the two electron signals. We now discuss how rhic were 
simulated. 



Derivation of the Result 



In central collisions of heavy projectile-target combinations of nuclei, for 
Ei a b > 1 GeV/u, we expect the formation of thermally equilibrated — or 
nearly equilibrated — matter (i. e. a "fireball") with subsequent hadroniza- 
tion; for beam energies below 1 GeV/u, no fireball should form, since in that 
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case heavy-ion collisions are adequately described by the participant specta- 
tor model, in which individual nucleons collide while the spectator nucleons 
remain cold and undetected. If E\ a \, is increased, the initial temperature of 
the fireball also increases, and then subsequently drops down as the fireball 
progressively expands. QCD should undergo a first order phase transition 
at a given temperature T c . For gluonic matter, lattice gauge theory suggests 
a transition temperature T c of about 207 MeV. For the temperature ranges 
involved, it is sufficient to consider QCD with just two massless quark fla- 
vors (u and d). In describing the time evolution of the strongly interacting 
matter, we take as the equations of motion relativistic hydro dynamics [|]. 

The picture of a hadronic fireball described by 1-D relativistic hydro- 
dynamics, with longitudinal expansion and subsequent pion freeze-out, has 
been experimentally verified by the CERN SPS central S + S collisions at 200 
GeV/c per nucleon||. This formulation requires equations of state (EOS), 
the notion of which implies the existence of an equilibrated initial uniform 
piece of nuclear matter. Since pionization of nuclear matter occurs already 
for T > 50 MeV, the EOS for hadronic matter in the relevant temperature 
range is chosen to be that of a (relativistic) pion gas plus a (non-relativistic) 
nucleon gas. The EOS for hadronic matter (pressure Phad vs. energy density 
phad) is given in parametric form by: 



where we denote = m n /T, fi n = n n +/T, etc. In each expression, the first 
sum indicates the contributions from 7r~, 7r° and 7r + , respectively (we use 




(1) 
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the fact that = —jj, n -, since ir + and 7r~ are mutually antiparticles, while 
/V) = since 7r° is its own antiparticle), and the second sum indicates the 
contributions from the nucleons (p = proton, n = neutron). For numerical 
purposes, we found it a good approximation for the EOS of the hadron phase 
to use Phad = 0.3 phad in the rhic temperature regime. 

For the quark-gluon plasma, the EOS is determined from a non-interac- 
ting relativistic gas of massless u and d quarks and gluons, along with the 
vacuum pressure (B constant), which takes into consideration color confine- 
ment as in a Bag model (this parameter is necessary for the phase transition 
to be of first order). The corresponding EOS is thus described by: 

P QGP ^f«|« 3 V ( — i + — i } + — T 4 - B (3) 

Q ir 2 Jo q ^ d I e x -^ + 1 e x+ ^ + 1 J 45 y ' 

PQ GP = %r dxx * E ( - 1 n + J t ) + ^ t4 + B ( 4 ) 
1 ^ 7T 2 Jo ^ . I e x -^ + 1 e x+ ^ + 1 J 15 v 7 

q=u,d 

where /i 9 = fJ> q /T, for q — u,d. The first and second terms in the sum are 
the contributions from q and g, respectively, while the last two terms are the 
contributions from the gluons and the (de-confined) vacuum, respectively. 
In the case of zero chemical potentials (zero net baryon number and electric 
charge), Eqs. (3) and (4) assume the simple form: 

Poor = ^7* -B (5) 

Pqgp = 37(^)^ 4 + B . (6) 
which, for T c = 207 MeV give B 1 / 4 = 288 MeV. 

In the general case of a net baryon number and/or net electric charge 
present, the equilibrium between the two phases is determined by the three 
equations in the three unknowns T c , fj, c , Pb (critical temperature, electric 
charge chemical potential, baryonic number chemical potential): 

7T+ ? Ppi H>ni 

Q = M p + K+-K-+\(M u -Mn)-\(M d -M i ) (7) 
B = Af P + Af n + l^u-Mu) + l(Md-m 
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where Q, B, Af p , Af n , Af u , Afu, • • • , etc are the respective electric charge, 
baryon number, proton, neutron, u-quark, anti-u-quark,. . . , etc number den- 
sities. The electric charge and baryon number, Q and B, are input parameters 
and at the same time the two conserved quantities in the thermal problem, 
so there are only two independent chemical potentials, namely pc and pb, 
in terms of which we can express those for each species: 

/^+ (= = He 

Hp = Hc + HB 

Hn = fJ>B (8) 

2 1 

Hu (= ~Hu) = ^Hc + ^Hb 

1 1 

Hd (= -m) = ~^Hc + ^Hb 

In the simultaneous solution of these equations, we find that the effect of a 
net baryon number on the value of T c is negligible, for the densities used in 
the rhic simulation. Thus the important energy densities, which are shown 
in Table 1, are fixed by the temperature alone. 

As the collision energy is increased, the first appearance of the QGP 
occurs when the energy density reaches p^id (0-22 GeV/fm 3 ), at T c = 207 
MeV. Because the phase transition is modeled to be first order, a latent 
energy density, ~ 4B, must then be supplied in order to liberate the quarks 
and reach a 100 % QGP. This happens when the energy density reaches Pqqp 
(3.81 GeV/fm 3 ). If the initial energy density, p i: is between 0.22 and 3.81 
GeV/fm 3 , the matter is in a mixture of the two phases, in which the initial 
proportion of QGP present, fc, satisfies p { = Pq GP x /• + Phld x (1 - fo)- 

Two parameters describe the collision previous to thermalization: the 
compression factor, k, and the stopping power, S. The volume V/j, of the 
initial fireball is taken to be the sum of the volumes of the individual nuclei 
divided by the compression factor: 

V fb = v n x (Ai + A 2 )/k, v n ~ 7.2 fm 3 /nucleon, (9) 
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where v n is the nuclear specific volume, and A\ and A2 are the mass numbers 
of the colliding nuclei. The old idea that the fireball can be infinitesimally 
small, based on the fact that the colliding nuclei are Lorentz contracted 'pan- 
cakes', is rejected; modern computations disclose that a finite time interval 
is required for equilibration and the compression factor k reaches a quasi- 
constant value of the order of 4.0 at the highest energies @. Thus we have 
adopted a sliding value of k from 1 to 4. The other major parameter is the 
stopping power, S, i. e. the average fraction of CMS kinetic energy a nucleon 
loses in the collision. For identical nuclei of baryon number A, the CMS 
thermal energy available for hydrodynamic expansion is 

E T = 2Ax (E N - m N ) x S, (10) 

where En is the CMS incident energy per nucleon and uin the nucleon mass. 
Experimentally, S ~ 0.7 only at high energies, while modest for small en- 
ergies. Thus we also made S a sliding linear function of En. These two 
parameters are set by the CERN SPS S + S rhic data|| to be about S ~ 0.5 
and k ~ 1.4 for incident CMS energy per nucleon En = 9.7 GeV. 

The last consideration is the choice of the relativistic hydrodynamic de- 
scription. Event shapes at high energy are not isotropic, but rather elon- 
gated along the beam axis|J so we choose a 1-D longitudinal simulation. A 
phenomenological parton Monte-Carlo simulation |7| gave results remarkably 
resembling those of a 1-D longitudinal hydrodynamic flow@]. 

Full numerical solutions including the EOS (l)-(4) and phase coexistence 
conditions (7)- (8) were studied. As mentioned above, effects of net baryon 
charge were found not to be significant, in which case the approximate EOS 
Phad = -Sphad for hadrons and for Eqs. (5)- (6) for QGP matter allow an 
analytical solution to the 1-D hydro flow: let r = \Jt 2 — z 2 be the proper 
time of a given fluid element (t,z being the spacetime coordinates from the 
origin of collision in the CMS); using r = 1 Fermi/c as the initial time 
{i.e. the estimated time for formation and equilibration of the fireball), for 
a high initial temperature the fluid element will start as QGP, cool down as 
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it expands reaching T c at r = Ti, condense to a full hadronic phase at r = r 2 
and subsequently cool below T c until freeze out. The temperature T and 
energy density p in the QGP, mixed and hadronic phases, respectively, are 
given by: 

QGP : T(r)/T(r ) = (ro/r) 1 ^ (p(r) - B)/(p(r ) - B) = (r /r) 4 / 3 
mixed: T(r) = T c , p ( r) = p QGP x f{r) + Phad x (1 - f{r)) (11) 
hadronic : T(r)/T(r 2 ) = (t 2 /t) - 33 , Phad (r) / Phad fo) = (r 2 /r) L3 

where /(r) = 1 - [1 - {ti/t)}sq GP /(sq G p - s had ) is the fraction of QGP in 
the mixed phase, s being the entropy density s = (p + P)/T. 

The simulation now proceeds as follows: at each beam energy, the initial 
energy density and temperature are calculated. Determination is made if the 
strongly interacting matter is in the hadronic phase (pions), mixed phase, or 
pure QGP phase. An estimate of the prompt electrons is then done, by com- 
puting the di-lepton rate for one time step (units in Fermi/c) over the initial 
fireball volume of a pion gas (the prompt electrons come from the initial fire- 
ball formation and do not know the existence of any subsequent QGP). After 
this, the hydro-expansion begins: the edges are advanced one time unit using 
the hydro-solution and the fluid is divided into zones. The new temperature 
is recomputed for the middle of each zone and the appropriate di-lepton rate 
determined, taking into account that the zone may be in a new phase. The 
electron rate in the QGP is due to quark-quark and quark-gluon collisions [§] 
(Fig. 3) and in the hadronic phase is due to pion annihilation || and virtual 



bremsstrahlung from pion scattering [TO]. Here we did not include heavier 



mesons JTTJ in our analysis, so that we may be slightly underestimating the 
lepton yield from the hypothetical hadronic fireball at high temperatures (i.e. 
above lGeV, that is, well above T c ). Also, finite temperature and density 



effects were not included, since they are small corrections [0 to the pion an 



nihilation amplitude in vacuum. Those zones which have an energy density 
less than half-nuclear experience freeze-out and do not produce any more 
particles (electrons and mesons). The approximate freeze-out temperature 
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is 154 MeV, with the approximate pion number density .135 fm~ 3 . By inte- 
grating over the invariant mass of the lepton pairs and then summing each 
zone contribution, the total production of electrons is ascertained for that 
time interval. We consider a fireball with mirror symmetry around z = and 
cylindrical symmetry about the beam axis. The time counter is advanced an 
additional Fermi/c and the hydro-expansion continues until all zones expe- 
rience freeze-out. The result, Fig. 1, for the number of produced electrons 
as a function of the beam energy is obtained. As mentioned earlier, the 1-D 
hydrodynamic model works very well with the CERN SPS S+S meson data, 
which constrains the model parameters k and S. For those values, Fig. 2 
shows the pion yield, corresponding to the number of pions in the fireball at 
the moment of freeze-out. 

There is a distinctive increase in the number of directly produced elec- 
trons, once the beam energy becomes high enough for the QGP phase to 
form. The physical reason for this manifold increase of directly produced 
electrons due to the existence of the QGP is the increased number of de- 
grees of freedom. This signal is independent of the particulars of the rhic 
simulation. Indeed, for the broad range of values of the phenomenological 
parameters considered, namely the stopping power S, the compression factor 
k and the effective gluon mass, the resulting direct electron signal showed 
no appreciable change. Thus, one might also expect an increased number of 
hard photons, although the analysis of that signal was not considered in this 
work. 
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Figure Captions 

1. The number of electrons per collision versus beam energy (GeV per 
nucleon), for Au vs. Au, after subtraction of the Dalitz noise. The 
CMS energy is twice the beam energy, for identical nuclei. 

2. The number of pions per collision versus beam energy (GeV per nu- 
cleon), for Au vs. Au. The CMS energy is twice the beam energy, for 
identical nuclei. 

3. The processes contributing to the production of dielectrons in the 
quark-gluon plasma. 
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Table 



e + e freeze-out (~ half nuclear energy density) 0.065 GeV/fm 3 

nuclear energy density 0.13 GeV/fm 3 

hadron energy density at the phase transition: p^ ad 0.22 GeV/fm 3 

QGP energy density at the phase transition: Pq GP 3.81 GeV/fm 3 



Table 1: Relevant energy densities in rhic. 
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Beam Energy (GeV/nucleon) 



